Effects of QCD phase transition on gravitational radiation 
from two-dimensional collapse and bounce of massive stars 
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We perform two-dimensional, magnetohydrodynamical core-collapse simulations of massive stars 
accompanying the QCD phase transition. We study how the phase transition affects the gravitational 
waveforms near the epoch of core-bounce. As for initial models, we change the strength of rotation 
and magnetic fields. Particularly, the degree of differential rotation in the iron core (Fe-core) is 
changed arametrically. As for the microphysics, we adopt a phenomenological equation of state 
above the nuclear density, including two parameters to change the hardness before the transition. 
We assume the first order phase transition, where the conversion of bulk nuclear matter to a chirally 
symmetric quark-gluon phase is described by the MfT bag model. Based on these computations, 
we find that the phase transition can make the maximum amplitudes larger up to ~ 10 percents 
than the ones without the phase transition. On the other hand, the maximum amplitudes become 
smaller up to ~ 10 percents owing to the phase transition, when the degree of the differential rotation 
becomes larger. We find that even extremely strong magnetic fields ~ 10 17 G in the protoneutron 
star do not affect these results. 

PACS numbers: Valid PACS appear here 

INTRODUCTION 

It was presented that quark matter may exist in the 
universe [l| . The existence is now supposed to be inside 
the core of neutron stars or inferred to be bare quark 
stars, where several observational signals have been sug- 
gested @, @, 0, @| . Much attention has also been paid 
to explore the relevant astrophysical phenomena in such 
sites d, 0, S, H, E3] ■ Related to the newer observation 
of PSR J0751+1807, skyrmion stars have been proposed 
to explain the large mass (~ 2.1 Mq) of a compact ob- 
ject Ei. 

It has been presented that the quark matter might ap- 
pear during supernova explosions [H, fl3l) or the tran- 
sition of a neutron star to a quark star [14[. Current 
supernova studies demonstrate that the stellar collapse 
of stars below ~ 25M in the main sequence stage leads 
to the formation of neutron stars, while in the case of 
more massive stars, to the formation of the black hole 
(l5| . In the latter case, quark matters might appear be- 
cause the central density could exceed the density of the 
QCD phase transition. It is worth mentioning here that 
supernova models with the phase transition are impor- 
tant because the vast release of the gravitational energy 
at the transition could be responsible for some classes 
of long-duration gamma-ray bursts [Io| . In this connec- 
tion, another energy source has been recently proposed 



core-collapse supernovae with the phase transition. We 
might study such microphysical phenomena in the core 
of the star directly from the detections of the gravita- 
tional waves. Currently, the gravitational astronomy is 
now becoming a reality. In fact, the ground-based laser 



interferometers such as TAMA300 18|, ll9( and the first 
LIGO [20, Hll are beginning to take data at sensitivi- 
ties where astrophysical events are predicted. For the 
detectors including GEO600 and VIRGO, core-collapse 
supernovae especially in our Galaxy, have been supposed 
to be the most plausible sources of gravitational waves 
(see reviews, for example, [22l l23j|). 

So far, there has been extensive work devoted to study- 
ing gravitational radiation in the context of rotational 



24, 
32, 



25, IgJ, [271 |28l l2f| [30, |3l( and magnetorotational 



33j core collapse supernovae. On the other hand, 



as quark nova [16l.ll7| 



It should be noted that the uncertainty of the equation 
of state (EOS) is always a big problem in the research of 



there are a very few simulations concerning the effects of 
the phase transition on the gravitational signals. Lin and 
his collaborators recently presented a three-dimensional 
simulation with the use of the polytropic equation of state 
in baryonic phase and MIT bag model in quark phase. 
They investigated how the delayed collapse of a rapidly 
rotating neutron star induced by the phase transition 
long after its formation, could produce the gravitational 
waveforms [34J . They demonstrated that the waveforms 
can be characterized by the damping timescale of the core 
and the induced core oscillation frequency. Furthermore, 
they pointed that the energy release in the form of the 
gravitational radiation owing to the transition could be 
less than ~ 10% of the gravitational binding energy. 
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In the present paper, we also pay attention to the ef- 
fects of the phase transition on the gravitational wave- 
forms, not at the epoch considered above, but at the mo- 
ment of core-bounce during the gravitational collapse of 
the massive stars. For this purpose, we perform the two- 
dimensional magnetohydrodynamic (MHD) simulations 
of the supernova cores accompanying the phase transi- 
tion. To treat the QCD phase transition, we employ 
a MIT bag model. We construct the initial models by 
changing the strength of rotation and magnetic fields, 
and the degree of differential rotation parametrically. We 
also vary the hardness of EOS in the baryonic phase in 
a parametric manner. In this paper, we investigate the 
relation between the phase transition and the change in 
the gravitational signals systematically. 

This paper is organized as follows. In sec. II we outline 
the numerical methods, input physics, the initial models, 
and the numerical methods for the gravitational wave- 
forms. In Sec. Ill, we present numerical results. Sec. IV 
is devoted to the conclusion. 



NUMERICAL METHOD AND INITIAL MODELS 
Numerical method 

The numerical method for MHD computations em- 
ployed in this paper is based on the ZEUS-2D code [IBj], 
(see [32| for details of the application to the core-collapse 
simulations). In the following equations, geometric units 
are used, G = c = h = 1. The basic MHD equations are 
written as follows, 



where 



Dp 
Dt 



-pV ■ v = 0, 



(1) 



Dv 1 
P —=-VP-pV^ cS +—(VxB)xB, (2) 



4-7T 



— = V x x B), 




A$ = Airp. 



(3) 

(4) 
(5) 



Here, p,P,v,e, B, <£>, and D/Dt, are the density, pres- 
sure, velocity, internal energy density, magnetic field, 
gravitational potential, and the Lagrangian derivative, 
respectively. In addition to the previous version [3^ |. 
we newly take into account the general relativistic cor- 
rection to the Newtonian gravity because it affects the 
waveform of the gravitational wave at the moment of the 
phase transition and the subsequent contraction of the 
core. The effective potential $ ff in Eq. © includes the 
correction (36|, which is defined to be 

$off = $ + 5$tov, 



(5$ TOV = 0> xov - $ s . 



Here, ^tov is the gravitational potential which is con- 
structed using the Tolman-Oppenheimer-Volkoff equa- 
tion. $s is the spherical Newtonian gravitational poten- 
tial. Albeit simple, this method has been demonstrated 
to reproduce well the results of the general relativistic 
calculations [3a| . 



Equation of State 

To describe the QCD phase transition in the super- 
nova cores, we follow the method adopted in Ref. [12| . 
The transition to a hadronic phase is modeled by adding 
the MIT bag constant B to the energy density. This 
phenomenological term lowers the pressure of the quark 
phase, where a transition to a confined phase occurs. 

The EOS in the quark phase is written as : 



i/3 r 



P= - 




(6) 



(7) 



(8) 



Here, p and pb are the energy density and the mean 
baryon chemical potential, respectively. We assume that 
Pb is three times that of the quark chemical potential. 
The number of quark flavors Nf — 3, and rib is the mean 
baryon number density which is assumed to be the one 
third of the quark number density. Using (O and ([7]), we 
get : 

p = 3P + 4B. 

We assume that the first order phase transition occurs 
during the collapse beyond some critical density, which 
has been supported by a recent study about the quark- 
gluon plasma [37|. We set the QCD coupling constant 
a s = in most models, which indicates that the coexis- 
tence phase between the baryon and quark is the widest 
in the density region. This choice helps us evaluate the 
effects of the phase transition maximally [l2j ]. In the 
baryon phase, we adopt the phenomenological EOS of 
BCK |38| that includes two parameters of the incom- 
pressibility K and the adiabatic index T to express the 
degree of the hardness of matter above po. The EOS of 
BCK is parametrized as follows, 



Kn 



(9) 
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Here the chemical potential in the baryon phase is given 



P 



M = 



rib 



(10) 



here no is the saturation number density (no = 
0.17 fm -3 ), and po = m u no (m u is the atomic mass unit). 
Two parameters of K and T are shown in Table |TJ We 
note that recent measurements give the constraints on 
the incompressibility : 210 < K < 270 MeV [2i3, 13 El • 
Therefore the adopted value of K = 220 MeV is consis- 
tent with this measurement. We furthermore take into 
account the parameter of K = 375 MeV, which would 
correspond to an extreme value. We add the pressures 
of electrons and photons to the EOS of BCK [12J. Since 
we focus on the behaviour of the EOS above n , we use 
the phenomenological EOS of Ref. 42j below no m ah 
models. 

As for the criterion of the phase transition, we im- 
pose the Gibbs condition with respect to the pressure 
and chemical potential to bridge the baryon and quark 
phase. We assume that the transition starts at p\ = 
5 x 10 14 g cm -3 for the baryon phase in most mod- 



whose value lies between those adopted by Gentile 
12j. Then, the end point of the transition (pi) 



els 
et al 

and B are determined analytically from the condition of 
P = Pi, p = fii, and a s = or 0.25 for the quark phase : 
pi and/or B(p\,P\, p\,a s ). Thus, we obtain p2, B, and 
Pb for two parameters of K and T as shown in Table [T] 
The bag constants used in this paper are consistent with 
the values used to investigate the structure of hybrid and 
quark stars [43j, |44 1 . 

While we get a reasonable value of pb ~ 950 MeV, 
pi and p2 are rather low compared to the suggestion 
from the QCD theory based on the lattice QCD cal- 
culations. It is predicted that the transition occurs at 
Pb > 1000 MeV and p2 > 10 15 g cm" 3 corresponding to 
the density of a neutron star [45j . We note that the pres- 
sure is not constant during the phase transition, since the 
electron and the photon pressures are included. 



Initial models 

We adopt the presupernova model of 13 Mq that has 
the iron core (Fe-core) of 1.2 Mq, in most models. The 
calculational area extends to 4000 km from the center 
with the mass of 1.4 M Q included. FIG.Q]shows the dif- 
ference of the pressure between the original presupernova 
model and ours whose EOS is described in Sec. II B. We 
can see that our EOS is phenomenological but consistent 
with original one. To examine the mass dependence of 
presupernova models, we adopt an another presupernova 
model whose mass is 40 Mq with the Fe-core of 1.9 Mq. 
In this model, the calculational area of 4000 km corre- 
sponds to 2.4 Mq. 




0.4 0.6 0.E 
mass [solor mass] 

FIG. 1: The difference of pressure between the original presu- 
pernova model (PSN) and ours (Our PSN), where the region 
of 1.4 Mq from the center is shown. 



Since the effects of the magnetic field and the angular 
momentum distribution on the presupernova star are un- 
certain (4t3 | , we make precollapse models from the nonro- 
tating progenitor model by adding the angular momen- 
tum and the magnetic field according to the following 
prescription. For the initial angular velocity distribution, 
we adopt the shell — type rotational law [461 ]: 



fi(r) = Q a 



where f2(r) is the angular velocity and r is the radius. 
Both fio and Ro are model constants that prescribe the 
rotational law. We regard the model with Rq = 10 3 km 
as uniformly rotating, since the radius of Fe core is ~ 10 3 
km. In most initial models, the initial magnetic field is 
assumed to be constant, Bo, which is poloidal and par- 
allel to the rotational axis in the computational domain. 
For the only two models in Table [Til we assume that the 
initial magnetic field is purely toroidal: 



B+ir) = Bo x 



R l 



where B^ (r) is toroidal component of the magnetic fields 
and Bo is a constant. The dominance of the toroidal 
fields to the poloidal ones ispredicted by the recent stellar 
evolution calculation 15|, |47J ■ As is discussed in Sec. Ill, 
effects of the magnetic field distribution on the transition 
are found to be small. 

We perform core-collapse simulations of 27 models 
which are given in Table [HJ The characters in the left 
hand side for each column indicate the initial condition 
concerning the magnetic fields and rotations: S (spheri- 
cal), u (weak uniform rotation and weak magnetic field), 
U (strong uniform rotation and strong magnetic field), 
R (strong differential rotation), d (strong differential ro- 
tation and strong magnetic filed), B (strong differential 
rotation and strong magnetic field of the shell— type law), 
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TABLE I: Physical quantities for EOS in hadron/quark phase equilibrium. The values, a„, pi, P2, B 1 ^, and pb are the QCD 
coupling constant, the critical density of the baryon phase, the density of the corresponding quark phase, the bag constant, and 
the mean baryon chemical potential, respectively. K and F are the incompressibility and the adiabatic index, which change 
the hardness of the baryon phase. 



Model 


a B 


Pl (10 14 g cm" 3 ) 


p 2 (10 14 g cm" 3 ) 


B 1 /4(MeV) 


MMeV) 


K(MeV) 


r 


-Mh 





5 


7.34 


164.8 


967 


375 


3 


-Mm 





5 


7.01 


163.8 


949 


220 


3 


-Ms 





5 


6.93 


163.8 


949 


220 


2.5 


-Mma 


0.25 


5 


6.00 


157.3 


953 


220 


3 


-Mm6 





6 


7.28 


163.8 


964 


220 


3 



D (very strong differential rotation and strong magnetic 
field). If there is "40" in the model name, it means that 
the mass of the presupernava model is 40 Mq . The char- 
acters after the hyphen indicate the adopted EOSs. The 
inclusion of the QCD phase transition is indicated by 'M': 
B (BCK without the phase transition), and M (MIT bag 
model with the phase transition). The right hand side 
characters mean the hardness for the EOSs in the baryon 
phase or another choice in the transition parameters. 
This hardness is expressed by the adiabatic index F and 
the incompressibility K: h (hard ; V = 3, K = 375 MeV), 
m (medium ; V = 3, K = 220 MeV), and s (soft ; T = 
2.5, K = 220 MeV). Two models u-Mma and d-Mma in- 
dicate that a — 0.25. The model u-Mm6 means that the 
critical density in baryon phase is p\ — 6 x 10 14 g cm" 3 . 
We do not calculate the case which corresponds to the 
model d-Mm, because the rotation of the core is so fast 
that the maximum density does not reach the critical 
density. 



with p being the matter density. The amplitude (jTTJ) is 
then transformed to the spherical coordinate as 



El 



h TT = Hi? = 



TT 



15 \ . 2 A 20 

— sin a 

TT I R 



(13) 



where a is the angle between the symmetry axis of the 
source and the direction to the observer and A E q is the 
second derivative of the radiative quadrupole , 



rl 2 

A E2 _ a M E2 



Af 2 consists of the following two terms: 



A E2 _ aE2 , A E2 

A 20 — ^20, quad ' ^20, Mag- 



(14) 



(15) 



Here, -A^j 2 j is the contribution from the matter. The 
magnetic component in Eq. (|15p is decomposed into two 
terms: 



Gravitational wave formulae from the rotating 
magnetized stellar cores 

To compute the gravitational waveforms from the ro- 
tating and magnetized stellar core, we follow the method 
of the quadrupole formula derived in [32|]. We describe 
the gravitational wave amplitude (GWA) hij that is cal- 
culated by 



2 d 2 TT , 
Rdt 2 ij [ 



R), 



(11) 



where subscripts i and j take over x, y, and z. t is the 
proper time and R is the distance from the source to 
an observer, respectively. The superscript TT indicates 
the transverse traceless part of the metric. is the 

reduced quadrupole moment, defined as 



rTT 



where is the total energy density including the contri- 
butions from the magnetic field, 

B 2 



(12) 



.4 



E2 

20. Ma 



A E2 



A E2 



(16) 



where A® 2 ■ 
and A^ 



kB is the contribution from the j x B part 
is that from the time derivatives of the en- 
ergy density of the electromagnetic field. In consequence, 
Eq. (fl3|) is composed of three terms : 



, TT 



,TT 
"quad 



n jxB 



(17) 



(see Ref. [32J] for details). In the following, we assume 
that the distance from the observer (R) is 10 kpc in the 
direction of the equator (a = ir/2). 



NUMERICAL RESULTS 

We summarize the physical quantities obtained from 
the numerical simulations in Table IIIIl As a guide to 
see the effects of the phase transition on the maximum 
amplitudes, we prepare the following quantity, 



A\h 



TT I 



I7) TT I 
\ ll M |max 



l/) TT l 

I ' L B I max 



(18) 



B I max 
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TABLE II: Model parameters. This table is divided into three groups : spherical (S- models), uniform rotation (u, 40u, U- 
models) and differential rotation groups (d, R, 40d, B, D- models). Note that the models with second capitals of "B" or "M" 
correspond to those without or with phase transition, respectively. 



Model 




Pi 


K 


r 


Ro 


T/\W\in\ 


E m /\W\im 


n 


Bo 






^ia14 „ ...... — 3\ 

ylu g cm ) 


(JVieV ) 




(km) 


(%) 


(%) 


( s ) 


(G) 


S-Bm 


- 


5 


220 


3 


- 


- 


- 


- 


- 


S-Mm 





5 


220 


3 


- 


- 


- 


- 


- 


u-Bh 


- 


5 


375 


3 


10 3 


0.1 


10~ 3 


2.3 


6.7 x 10 IU 


u-Mh 





5 


375 


3 


10 3 


0.1 


lO" 3 


2.3 


6.7 x 10 10 


u-Bm 


- 


5 


220 


3 


10 d 


0.1 


10" 3 


2.3 


6.7 x 10 1U 


u-Mma 


0.25 


5 


220 


3 


10 3 


0.1 


IO -3 


2.3 


6.7 x 10 10 


u-Mm6 





6 


220 


3 


10 3 


0.1 


IO -3 


2.3 


6.7 x 10 10 


u-Mm 





5 


220 


3 


10 3 


0.1 


IO" 3 


2.3 


6.7 x 10 10 


u-Bs 


- 


5 


220 


2.5 


10 d 


0.1 


10 a 


2.3 


6.7 x 10 IU 


u-Ms 





5 


220 


2.5 


10 3 


0.1 


IO -3 


2.3 


6.7 x 10 10 


40u-Bm 


- 


5 


220 


3 


10" 


0.1 


10 3 


2.1 


9.2 x 10 IU 


40u-Mm 





5 


220 


3 


10 3 


0.1 


IO -3 


2.1 


9.2 x 10 10 


U-Bm 


- 


5 


220 


3 


10 


0.5 


i n— 1 

10 


5.3 


6.7 x 10 11 


U-Mm 





5 


220 


3 


10 3 


0.5 


10 _1 


5.3 


6.7 x 10 11 


d-Bm 


- 


5 


220 


3 


10* 


0.5 


10" 1 


58.8 


6.7 x 10 11 


d-Mma 


0.25 


5 


220 


3 


10 2 


0.5 


io- 1 


58.8 


6.7 x 10 11 


d-Mm 





5 


220 


3 


10 2 


0.5 


lO" 1 


58.8 


6.7 x 10 11 


R-Bm 




5 


220 


3 


10 2 


0.5 




58.8 




R-Mm 





5 


220 


3 


10 2 


0.5 




58.8 




40d-Bm 




5 


220 


3 


10* 


0.5 


IO" 1 


79.4 


9.2 x 10 11 


40d-Mm 





5 


220 


3 


10 2 


0.5 


io- 1 


79.4 


9.2 x 10 11 


B-Bm a 




5 


220 


3 


10 2 


0.5 


lO" 1 


58.8 


1.3 x 10 14 


B-Mm 





5 


220 


3 


10 2 


0.5 


io- 1 


58.8 


1.3 x 10 14 


D-Bm 




5 


220 


3 


50 


0.5 


IO" 1 


183.9 


6.7 x 10 11 


D-Mm 





5 


220 


3 


50 


0.5 


io- 1 


183.9 


6.7 x 10 11 



"For two models of B-Bm and B-Mm, the shell-type distribution 
of the magnetic field B^ = Bq x R^/(r 2 + Rq) is assumed. 



where l/i^/lmax and |/ig T | ma x are the absolute value of 
gravitational wave amplitudes (GWAs) with and without 
the phase transition, respectively. 

From the table, we can see that the phase transition 
makes the maximum GWAs larger by a few up to ~ 10 
percents for the uniformly rotating models. On the other 
hand, the phase transition lowers the maximum GWAs 
for the differentially rotating models. If the rotational 
strength T/|W| becomes very large, the effects of phase 
transition on gravitational waves become small. We ex- 
plain these features in the following subsections, where 
the changes in GWA from models to models are exam- 
ined. Furthermore, we perform the Fourier transforma- 
tions for all models to show more comprehensive and sys- 
tematic analysis of the gravitational waveforms. 



formly rotating models as seen in Table HTT1 and FIGs. 

As the initial rotational strength T/|W|ini increases, 
the maximum GWA, |/i TT | m ax becomes large. However 
the effect of the phase transition on the maximum GWA 
decreases as seen from A|/i TT | max . In the right panels 
of FIGs. [2]|4l we can see that the frequency region with 
the phase transition shifts to the higher one as a whole, 
because the increased density during the collapse tends 
to the higher frequency [3J|. As for the effects of the 
equation of state, we find that the larger gravitational 
wave is radiated for the soft EOS, regardless of the phase 
transition. Comparing the models of hard, medium and 
soft EOS in Table ITITl this tendency is recognized to be 
independent on whether the softness originates from T or 
K. 



Effects of the phase transition and rotation 

We will first show the effects of the phase transition 
on GWA for the uniformly rotating models. We find 
that GWAs are larger from a few percents to about ten 
percents by considering the phase transition for all uni- 



To understand these results, we give an order-of- 
magnitude estimation of the quadrupolc formula 
The gravitational amplitude is proportional to the second 
time derivative of the radiative quadrupole (see Eqs. (I13p 
and HI])). The second time derivative in Eq. (fT4")) could 
be replaced by the reciprocal of the square of the time, 
if the time scale is very short. The typical dynamical 



6 



TABLE III: Physical quantities for the models with and without the phase transition, tb, ifln and p max are the bounce time, 
the final time of the calculation, and the maximum density. T/|W|fi n and i5 m /|W]a n are the final values of T/I1FI and i5 m /|W|, 
respectively. |/i TT | ma x is the maximum GWA, and A|/i TT j max is defined by Eq. (|18p . Note that models with second capitals 
of "B", "M" corresponds to without and with phase transition, respectively. fi st is the first peak frequency in the Fourier 
transformation of GWA. 



Model 


h 




Pmax 


T/\W\en 


£m/|!U| fln 


\h TT \ 


A\h TT \ 


fist 




(ms) 


(ms) 


(10 14 g/cm 3 ) 


(%) 


(%) 


(lO" 20 ) 


(%) 


(Hz) 


S-Bm 


75.5 


126 


6.95 












S-Mm 


75.5 


146 


8.84 












u-Bh 


76.3 


82.6 


5.37 


1.24 


2.12 xlO 


0.28 




727 


u-Mh 


7fi 3 


89 7 




1.24 


9 xin -3 


n 9Q 


iq a 


748 


u-Bm 


76.2 


82.7 


6.90 


1.29 


2.07 xlO -13 


0.31 




652 


u-Mma 


76.3 


87.2 


7.19 


1.31 


2.24 xlO" 3 


0.32 


+3.2 


652 


u-Mm6 


76.2 


83.4 


8.05 


1.27 


2.24 xl0~ 3 


0.33 


+6.4 


838 


u-Mm 


76.3 


87.0 


8.60 


1.33 


4.27 xl0~ 3 


0.34 


+9.6 


875 


u-Bs 


76.3 


83.3 


8.12 


1.29 


2.25 xl0~ 3 


0.34 




587 


u-Ms 


76.2 


86.6 


8.79 


1.36 


3.28 xl0~ 3 


0.38 


+11.8 


931 


40u-Bm 


70.6 


76.1 


7.02 


1.29 


4.15 xW~' ! 


0.26 




635(889) 


40u-Mm 


70.7 


79.4 


8.64 


1.39 


5.42 xlO -2 


0.28 


+7.6 


762 


U-Bm 


79.3 


86.0 


6.35 


5.05 


5.32 xl0~ 2 


1.25 




859 


U-Mm 


79.4 


93.2 


8.03 


4.97 


8.25 xlO -2 


1.27 


+1.6 


872 


d-Bm 


83.8 


90.9 


5.64 


9.41 


1.16 xl0 _i 


3.03 




792 


d-Mma 


83.8 


91.0 


6.27 


9.52 


1.15 xlO -1 


2.79 


-7.9 


743 


d-Mm 


83.8 


91.0 


7.31 


9.40 


1.14 xlO -1 


2.72 


-10.2 


659 


R-Bm 


83.8 


90.9 


5.67 


9.53 




3.03 




791 


R-Mm 


83.8 


91.0 


7.40 


9.51 




2.71 


-10.6 


658 


40d-Bm 


84.1 


89.9 


5.81 


12.8 


1.33 xlO" 1 


2.98 




562 
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time scale of the collapsing core is idyn ~ 1/ \/Gp. As 
the result, the amplitude is roughly proportional to the 
core density p c multiplied by the radiative quadrupole 
moment: 

h TT cx Pc Mf 2 . (19) 
The leftsides of Figs. OH show the time dependence of 



the maximum density. The maximum density of the 
model with the phase transition (for example the model 
of u-Mm) is always larger than that without the tran- 
sition (for example the model of u-Bm). In the figures 
of Pmax and their corresponding GWA figures, there are 
high frequency spikes. However, these are numerically 
ones, because quark matter areas (~ 10 km) are much 



2.0x10"' 



1.0x10 



21 



0.0x10" 

,-21 



-1.0x10"' 



-2.0x10 



21 



-3.0x10"' 



-4.0x10"' 









u-Bm 



0.074 0.075 0.076 0.077 0.078 0.079 0.08 
time [s] 



1.2x10' 




1 1.2 1.4 

frequency[kHz] 



FIG. 3: Same as FIG. but for u-Mm and u-Bm. 



2.0x10"' 



1.0x10"' 



0.0x10" 



-1.0x10"' 



-2.0x10"' 



-3.0x10"' 



-4.0x10"' 



A 


i i i 


2.0x1 0" 24 






1.6X10- 24 




I H / k/' 


f 1.2x10~ 24 
8.0x1 0~ 25 




u-Ms 

u-Bs 

i i i 


4.0x1 0" 25 



0.075 



0.076 



0.077 0.078 
time [s] 



0.079 



0.08 




FIG. 4: Same as FIG. [2] but for u-Ms and u-Bs. 



9.00x10 
8.00x10 
7.00x10 



14 



«~ 6.00x1 14 
-| 5.00x1 14 
*4.00x10 14 
cE3.00x10 14 
2.00x1 14 



1.00x10 
0.00x10 




0.074 



0.076 
time [s] 



-3.08x10"' 



-3.10x10"' 



-3.12x10"' 



■H, -3.14x10"' 



"J -3.16x10"' 



-3.18x10"' 



0.078 



-3.20x10"' 




0.072 



0.075 0.078 
time [s] 



0.081 



FIG. 5: Maximum density and radiative quadrupole moment of the model u-Mh and u-Bh. 



9.00x10 
8.00x10 
7.00x10 



14 



«~ 6.00x1 14 

-y 5.00x1 o 14 

CD 

*4.00x10 14 

J3.00x10 14 

2.00x1 14 
1.00x10 14 
0.00x10 




0.074 



0.076 
time [s] 



0.078 



-3.08x1 0" J 
-3.1 0x1 0" 3 
w -3.12x10" 3 
^ -3.14x10" 3 

O 
OJ 

^ -3.1 6x1 0" 3 
-3.18x10" 3 
-3.20x1 0" 3 



I I 
u-Mm 


i i 


_u-Bm 


/ \ 

/ V 




\ 

\- 

1 1 


i i 



0.072 



0.075 0.078 
time [s] 



0.081 



FIG. 6: Same as FIG. [5] but for u-Mm and u-Bm. 



9.00x10 
8.00x10 
7.00x10 
<^ 6.00x10 
-y 5.00x10 
*4. 00x10 
<^3.00x10 14 
2.00x10 
1.00x10 
0.00x10 




0.074 



0.076 
time [s] 



-3.08x10"' 



-3.10x10"' 



-3.12x10"' 



S -3.14x10"' 



-3.16x10"' 



-3.18x10"' 



0.078 



-3.20x10"' 



u-Ms 

ii-Rc 


i i 


N 

v- 
\ 

i i 



0.072 



0.075 0.078 
time [s] 



0.081 



FIG. 7: Same as FIG. [5] but for u-Ms and u-Bs. 



1.5x10"' 
1.0x1 CP 



5.0x10"' 
0.0x10" 



-5.0x10"' 
-1.0x10" : 



-1.5x10"' 
-2.0x1 0" : 



-2.5x10"' 
-3.0x1 0" : 




0.069 0.07 0.071 0.072 0.073 0.074 
time [s] 




1.5 

frequency [kHz] 



-3.5x10 



FIG. 8: Same as FIG.[2]but for models of 40u-Mm and 40u-Bm. 

2.5x1 0' 2 " 




0.075 



0.0755 0.076 0.0765 0.077 0.0775 0.078 
time [s] 




u-Mm 
u-Mma 



1 1.2 1.4 1.6 1 

frequency [kHz] 




FIG. 9: Same as FIG. [J but for the models of u-Mm and u-Mma. 

1.2x10" Z " 



1.0x10" 20 






5.0x1 0" 21 


" J\ 




0.0x10° 






5.0x1 0" 21 






1.0x10" 20 






1.5x10" 20 






2.0x1 0" 20 






2.5x1 0" 20 






3.0x1 0" 20 


V R-Mm 


3.5x1 0" 20 


R-Bm 





0.081 



0.082 0.083 0.084 0.085 0.086 0.087 
time [s] 




1.2 1.4 1.6 1.1 
frequency[kHz] 



FIG. 10: Same as FIG. [5] but for the models of R-Mm and R-Bm. 



9 



smaller than overall calculational areas (about Fe-core 
size). More fine tuned calculations with many meshes 
will delete such spikes. The absolute values of Mf 2 are 
almost the same around the peak, if the initial rotational 
strengths are the same (see the right panels of Figures [5} 
[7]). Since the transition causes the increase in the density 
being similar to the effect of soft EOS, we get larger am- 
plitudes. |/i TT |max in u-Ms increases by 10% compared to 
u-Bs as seen in Table IIIII On the contrary, we find that 
the effect of the phase transition becomes small if the 
rotational strength becomes as large as T/|W|fi n ~ 5%, 
which leads to the suppression of the contraction due to 
the strong centrifugal forces (compare p max of u-Mm with 
that of U-Mm in Table [m|. 

For the Fourier transformations, the ranges of frequen- 
cies are distributed in 500-2000 Hz for all models. Clear 
peaks are not found compared to the case in Ref. [34j |. 
While we adopt realistic presupernova models as initial 
models, they use polytropic models of equilibrium neu- 
tron stars. If QCD phase transition is considered, both 
the soft EOS model (FIG. [4| and slowly rotating model 
(FIG. [3]) result in the shift to the high frequency sides in 
proportion to the increased density during the collapse 
(see fist, of the models with/without the transition in 
Table 

To check the effects on GWA due to the difference of 
presupernova models, we calculate the collapse of 40 M Q 
models. The result is given in Table III (compare model 
40u-Mm with u-Mm, or 40d-Mm with d-Mm). We see the 
difference of the maximum density at the bounce, which 
is ascribed to the size of the Fe-core. The slight differ- 
ences of the maximum amplitudes are found between the 
models with/without the phase transition. The difference 
of the presupernova models dose not make a qualitative 
difference in both uniformly rotating and differentially ro- 
tating models. In the Fourier transformation, the overall 
frequencies for 40u-Mm shift to the higher region com- 
pared to the case of 40u-Bm as seen in FIG. [5]except for 
around fi st , where the first peak is not clear. 

We change a value of coupling constant a s . As the 
width of the density during the phase transition (coexis- 
tence area) becomes narrower for a s 0, the difference 
of the amplitude by the presence of the phase transition 
should become small. In Table [TTTl the value \h TT \ max of 
u-Mma is between u-Mm and u-Bm. In the Fourier trans- 
formation, the change in f\ st is small, but the frequen- 
cies shift overall between those of u-Mm and u-Bm (see 
Table [III] and FIG. E). In differentially rotating model, 
GWA and the frequencies of d-Mma lie between those of 
d-Mm and d-Bm, too. 

The difference of critical density (pi) for the gravita- 
tional wave is shown in Table HTT1 The higher critical den- 
sity (u-Mm6) results in the narrower width of the density 
region during the transition. As a consequence, \h TT \ max 
and fi st of u-Mm6 are between the corresponding values 
of u-Mm and u-Bm. 



For the strong differential rotation, it is shown in Ta- 
ble IIIII and FIG. [10] that the maximum amplitudes be- 
come smaller for the models with the phase transition 
(see A|/i TT | max of the three models, from the bottom in 
table IIIip . To explain such aspects, we refer to the early 
research [33j . Their results show that very fast rotation 
with soft EOS models tend to suppress the centrifugal 
force element of GWA. Here very fast rotation is not lim- 
ited in differential rotation, but includes very fast uni- 
form rotation whose T/lVFl is one digit larger than our 
uniformly rotating models. Since our initial models are 
spherical models, we do not adopt such a large T/|W| 
for consistency with our initial models. Consequently, 
the strong differentially rotating models with the transi- 
tion in our calculations are the same as very fast rotation 
with soft EOS, which lowers 1st peak of GWA for our 
differentially rotating models shown in Table III. Corre- 
sponding to these, their frequencies with the transition 
shift to low (see their Fourier transformation, FIG. [T0l 
and fist.s of Table iHl]). The above estimate (|T9|) is useful 
for the interpretation of the results. However it is found 
to be not applicable for strong differential rotation. This 
is because the differential rotation acts against a matter 
infall to the center by the phase transition due to the 
stronger centrifugal forces, but simultaneously leads to 
the stronger accretion to the center due to the smaller 
angular momentum in the rather outer regions. Due to 
this very subtle competition, we can only know from the 
numerical results that the differential rotation makes the 
amplitudes lower up to ~ 10 percents at the moment 
of the phase transition as far as the parameters in the 
present calculations. 

Effects of the phase transition and magnetic fields 

We focus on the models with the strongest magnetic 
fields ("B-" models in Table HJ) to clarify the effects of 
the magnetic field on the gravitational wave. First, we 
compare each component h^J ad , hJ^ B and tij^ of h TT in 
Eq. (fT"7|) . Figure [IT] shows the waveforms originated from 
the mass quadrupole moment, j X B part, and the time 
derivative of the magnetic energy density p m . The most 
definitive component to GWA is the one originated from 
the mass quadrupole moment h^J ad . 

To see the influence of the magnetic field on the gravi- 
tational wave frequencies, we compare the zero magnetic 
field models (R-Bm, R-Mm) with strong magnetic field 
ones (B-Bm, B-Mm). It is understood from FIG [TU]fT2l 
that the influence of the magnetic field with the tran- 
sition on the gravitational wave can be neglected. The 
corresponding maximum GWAs have almost the same 
values as shown in Table ITTTl 

The components of GWA for B-Bm and B-Mm are 
summarized in Table [TVl It should be noted that sign of 
each component at the bounce is different in the strong 
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FIG. 12: Fourier transformation of waveforms for the strongly 
magnetized model B-Mm and unmagnetized model R-Mm, 
corresponding to GWAs of FIG. [TO] and FIG. [TT] . 



TABLE IV: Three components in Eq. (JT7J- Note that the 
listed values correspond to the time when the values of 
I^Jxs/^mmd \max is maximum. It means that the time is not 
same with the time of |/i TT | m ax in table ITTTl All amplitudes 
are given in units of 10 -20 . 
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Bm (R-Mm) due to this contribution from the magnetic 
fields. 

At the moment of the phase transition, the magnetic 
fields in the central regions become larger due to the com- 
pression, because the magnetic fields are frozen-in to the 
matter. As a result, the contribution to the amplitudes 
from the magnetic fields become larger for the model with 
the transition, but the change is found to be as small as 
1% (compare |/iJJ B //iq U T ad | in table HE!) 



FIG. 11: Waveforms of the strong magnetized models, B-Mm 
and B-Bm. The top figure is one originated from the mass 
quadrupole moment, the next figure is from j x B part and 
the bottom is from the time derivatives of the magnetic energy 
density p m . Note the scale differences of the vertical axis. 

magnetic field models. As the result of cancellation be- 
tween h^ ad and hj^ B , the ratio of |/iJJ B //iq U T ad | becomes 
less than 8 % regardless of the phase transition, which can 
be roughly estimated by the ratio of the magnetic to the 
matter energy density at bounce, namely, 

Kotake et al. (2004) [32j have found the same effect 
though they did not include the phase transition. In 
fact, we can see from Table HITl that GWAs of the model 
B-Bm (B-Mm) is smaller compared to other models R- 



CONCLUSION 

We have performed two-dimensional axisymmetric, 
magnetohydrodynamical simulations for supernova cores 
accompanying the QCD phase transition. To elucidate 
the implications of a phase transition against a super- 
nova, we investigated how the phase transition affects the 
gravitational waveforms near the epoch of core-bounce. 
As for the initial models, we changed parametrically the 
strength of the rotation and the magnetic fields. As for 
the microphysics, we adopted a phenomenological equa- 
tion of state above the nuclear matter density, including 
two parameters to change the hardness of the matter be- 
fore the transition. To treat the QCD phase transition, 
we employed a MIT bag model. Based on these computa- 
tions, we showed that the phase transition can make the 
maximum amplitudes larger up to ~ 10 percents than the 
ones without the phase transition. On the other hand, we 
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found that the phase transition makes the maximum am- 
plitudes smaller up to ~ 10 percents, when the iron core 
rotates strongly differentially. It was confirmed that the 
strong magnetic fields themselves decrease the maximum 
amplitudes by less than 8% regardless of the transition. 
Even the extremely strong magnetic fields ~ 10 17 G in 
the protoneutron star do not affect the above features. 

Finally, we give some discussions and speculations 
based on the obtained results. We assume rather low 
density for the onset of the phase transition in compar- 
ison with a recently predicted EOS by the lattice QCD 
computations [45| . Thus the results here could be some 
extreme cases predicted by the phenomenological MIT 
bag model, and the change in the amplitudes due to the 
transition could be interpreted as an upper bound. 

Although the other initial presupernova models may 
change our analysis to some extent, the qualitative re- 
sults obtained here will not be different so much. This 
is because the structure of the iron core is similar, while 
the initial mass of the helium core increases with the 
progenitor masses [48| . For the presupernova model of 
40Mq that has the iron core of 1.9M Q , we have found 
that GWA at the bounce is only about 3 percents larger 
than that of 13M Q . Since we concentrate on the forms 
of the gravitational wave just after the bounce, the dif- 
ference in the helium core mass should not be important 
at all. 

If we could calculate the hydrodynamical evolution 
long after the core bounce especially in some failed core- 
collapse supernova models, we would observe the phase- 
transition of the protoneutron stars. We also think 
interesting to investigate the possible effects of color- 
superconductivity recently proposed [491 ]. As a next step, 
we are now going to investigate how the gravitational 
waves will be originated from such events in the context 
of magnetorotational core-collapse. 
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